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Abstract 

We study the screening of static magnetic and electric fields in massless quantum 
electrodynamics (QED) and massless scalar electrodynamics (SQED) at tempera- 
ture T. Various exact relations for the static polarisation tensor are first reviewed 
and then verified perturbatively to fifth order (in the coupling) in QED and fourth 
order in SQED, using different resummation techniques. The magnetic and electric 
screening masses squared, as defined through the pole of the static propagators, are 
also calculated to fifth order in QED and fourth order in SQED, and their gauge- 
independence and renormalisation-group invariance is checked. Finally, we provide 
arguments for the vanishing of the magnetic mass to all orders in perturbation 
theory. 
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1 Introduction 



The simplest and probably best known manifestation of many-body effects in electro- 
magnetic plasmas is Debye screening : the usual Coulomb potential between two static 
charges in a vacuum transforms, in the presence of a plasma, into a Yukawa potential, 
V(r) ~ e~ mDr /r. The scale is called the electric (Debye) screening mass. For a 
plasma at high temperature T (much larger than the electron mass), m 2 D = e 2 T 2 /3 to 
lowest order in the coupling constant |L]. 

The above-mentioned relationship between the static potential and the screening 
mass is usually established within the approximation of linear response (see [[J and section 
2), whereby one calculates the potential between two arbitrarily weak external static 
sources q\ and q<i separated by r, 

f d 3 p e ij5V qiq 2 f°° dp pe ipr 

[r) = gig 2 J ( 27 r) 3 j9 2 + n L (0,p) = 2^1-oo2m p 2 + IT L (0, p) 

with lli(0,p) = — n o(0,p), n o(0,p) being the static (po = 0) electric polarisation tensor 
(our conventions and notations are summarized in section 2). The integral over p may 
be performed by closing the integration path in the upper half of the complex p plane. 
When the external charges are widely separated (r — > oo), the behaviour of the potential 
is dominated by the singularity of D 00 (0,p) = — l/(p 2 + U L (0,p)) which lies closest to the 
real axis. To leading order, this is a pole occuring for p 2 = U^(0,p -»■ 0) = -e 2 T 2 /3, 

(2) 

Ho being the one-loop polarization tensor. 

One of the main objectives of this paper is to study, within the definition ( |1.1| ), 
corrections to this leading (~ e 2 T 2 ) Debye mass, and in the process clarify several issues 
which arise. We shall verify in a perturbative calculation that the dominant singularity 
of D 00 remains a pole on the imaginary axis. This allows us to define mo as the solution 

of i 

m 2 D = -U 00 (0,p) \ p2= _ ml . (1.2) 

We shall perform the computation in eq. (|1.2j ) up to order e 5 in usual (spinor) electro- 
dynamics (QED), and up to order e 4 in scalar electrodynamics (SQED). In both cases, 
this corresponds to the first two non-trivial corrections above the "hard thermal loop" 
approximation |4], ||, and is an extension of previous lower order calculations. As we shall 
see, the Debye mass obtained in this way is renormalisation group invariant to the re- 
quired order. General arguments also indicate the gauge-independence of propagator 
poles and this too will be demonstrated in our computations. At this point it is worth not- 
ing that, though we have introduced the screening mass via the potential above, the same 
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mass Trip controls the exponential decay of other interesting gauge-invariant correlators 
such as that of static electric fields (Ei(x)Ej(y)). 

Solving eq. ( |1.2| ) requires the knowledge of il o(0,p) for p ~ m ~ eT. Thus we are 
led to study first the simpler object n o(0,0), up to the order of interest. In massless 
QED, this quantity has been already computed to order e 5 , but only indirectly, by using 
the relation [0, §] 

n 00 (o, P ^o) = -e 2 |4r, (1.3) 

(where P is the plasma pressure and /i e is the chemical potential for the charged particles), 
together with an old result for the pressure ||. Notice that, because of the explicit factor 
of e 2 on the right-hand-side (r.h.s.) of eq. (|1.3| ), the result for the pressure is only required 
to third order (but for finite chemical potential). Still, since the derivation of eq. ( |1.3| ) is a 
formal one (see section 2.2), it is of interest to verify the identity explicitly in some cases 
to show that it is not invalidated, for example, by potential problems such as infrared 
divergences. This, we shall do in the case of QED: we shall compute the left-hand-side 
(l.h.s.) of ( |1.3|) to fifth order (and at zero /i e ), and we shall check that it agrees with 
known results for the r.h.s. In the case of SQED, Iloo(0, 0) has been previously computed 
up to order e 3 |9j (see also Ref. |1(J); we will present in this paper the order e 4 correction. 

Another objective here is to look into the screening of static magnetic fields. It is 
known that such fields are not screened in ordinary plasmas. Indeed, one may rely on 
Ward identities, together with the (exact) Dyson- Schwinger equation to show that 

11^(0,^0) = 0, (1.4) 

(see Sect. 2.2 for more details). However, in order to guarantee the absence of pertur- 
bative singularities beyond that at p 2 = in the correlator of magnetic fields, one needs 
the stronger result Uij(0,p — > 0) = 0(p 2 ). This will be verified explicitly up to fifth 
(respectively, fourth) order in perturbation theory for QED (respectively, SQED). The 
perturbative arguments can be extended to an all-order proof, to be detailed in Sect. 3.1 
in the case of QED. 

As in all higher order perturbative calculations at nonzero temperature, a systematic 
determination of requires a resummation of the large collective plasma effects. For the 
calculation of dynamical quantities, the resummation involves the procedure developed 
by Braaten and Pisarski [|J whereby one uses effective vertices and propagators obtained 
by dressing the bare quantities with "hard thermal loops" . The latter are the dominant 
parts of one-loop amplitudes and express the effects of Landau damping, Debye screening 
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and collective oscillations 0, |5|, IT|. However for the computation of static (zero external 
energy) Green's functions, like those to be discussed here, the calculations are most natural 
and convenient in the imaginary time formalism, without analytical continuation to real 
time. Then the general programme of Refs. |4[] reduces to the well known Gell-Mann- 
Bruekner resummation studied many years ago in the nonrelativistic context and 



subsequently extended to the relativistic regime g, [13], [14 . 

For static calculations in the imaginary time formalism (which we use exclusively 
in this paper), the resummation concerns only the internal lines with zero Matsubara 
frequencies (also referred as static lines), and consists in dressing these lines with the 
corresponding screening thermal masses. No vertex resummation will be needed: in 
QED, there are no vertex in which all lines are soft; in SQED, there are no hard thermal 
loops beyond the two-point functions 0, |TD[. The non-static internal propagators need 
not be resummed since the corresponding Matsubara frequencies ensure an infrared cut- 
off of the order of T, relative to which all the thermal corrections are perturbative. This 
relative simplicity of resummation for static quantities allows different approaches to 
the detailed calculations. At low orders one can in general perform the resummation of 
diagrams by "inspection" ||13|| . In static QED calculations, a "resummation by inspection" 
(let us call this method (a)) is even feasible at very high orders fl5| because the fermion 
lines are always hard in imaginary time and so do not require dressing. On the other 
hand, in theories with self-interacting bosonic fields more efficient and practical methods 
of higher order calculations are required which perform the equivalent resummation, of 
which there are at least three: (b) truncation of the full skeleton expansion^, [L4j]; (c) the 
use of rearranged lagrangians incoporating the screening masses |L6], £|, [17], [18]] and, (d) the 
use of dimensionally reduced effective lagrangians obtained by systematically integrating 



out of the heavy modes [2(| [H]] (and references therein). In this paper we will employ 
resummation by inspection, method (a), for the QED calculations and the truncation of 
Schwinger-Dyson equations (method (b) above) for SQED. By way of comparison, we 
also discuss the calculation in SQED using the effective lagrangian, method (d). Method 
(c)will not be used in this paper but a recent discussion at high orders may be found in 
Ref.|]|. 

As a result of the resummation, the perturbative expansion for the electric mass 
involves odd powers of the coupling strength, that is, it is not analytic with respect to e 2 . 
This becomes apparent in spinor QED only at the order e 5 , but is already manifest at the 
order e 3 in scalar QED||, as well as in QCD0. Odd powers of the coupling occur in the 
perturbative expansion since, after dressing the soft propagators with the corresponding 
screening masses, the relevant expansion parameter in the infrared is e 2 (T/m) ~ e, rather 
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then e 2 . From the point of view of its infrared behaviour, SQED is more interesting than 
spinor QED since it involves interacting bosonic fields. As a consequence, the scalar theory 
bears more ressemblence to QCD, by exhibiting some of the non-trivial IR structure of this 
latter theory, but in the (technically less involved) context of an abelian gauge structure. 
In this respect, SQED serves as a toy model for investigating resummations techniques 
to be eventually applied in high-temperature QCD (see, e.g., Refs. for such recent 
applications) . 

The plan for the rest of the paper is as follows. In Sect. 2 we display our notation 



and conventions, and derive the relations (|L1 ), (|1.3|) and (|1.4j) in an exact, but formal 



manner. Although most of the material in this section may be found scattered in the 
literature, we have collected it in one place to keep the discussion self-contained. In Sect. 
3 we study massless QED in perturbation theory. The highlights in this section are the 
direct calculation of m 2 D to order e 5 T 2 , discussion of its gauge- invariance, and an all-order 
proof for the vanishing of the static magnetic screening mass. Section 4 is devoted to 
scalar QED. The notable results obtained here are 11^(0, p — > 0) and m 2 D to order e 4 T 2 , 
and the vanishing of the magnetic screening mass to the same order. In Sect. 4.6, we 
re-discuss the results for SQED from the point of view of the effective three-dimensional 
theory for static fields. This sheds a new light on the resummation, and helps keeping 
track of the various terms in the diagramatic expansion. We conclude in Sect. 5 with a 
summary of results and some discussion. Some technical details omitted from Sect. 4 are 
collected in the Appendices. 



2 Notation and General Results 
2.1 Conventions 

We summarise here our conventions and notation. Unless otherwise stated, all calcula- 
tions from Sect. 3 onwards will be for the massless theories at zero chemical potential. 
Ultraviolet divergences are regulated by dimensional continuation (4 — > D = 4 — 2e) and 
renormalisation is via minimal subtraction. We employ the imaginary time formalism and 
denote the four-momenta by capitals, = (go,q), qo = iu n — innT, with n even (odd) 
for bosonic (fermionic) fields. The scalar product is defined with a Minkowski metric, 
so that Q 2 = g 2 — q 2 . The measure of loop integrals will be denoted by the following 
condensed notation: 
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f[dQ) = T £ / (dq) , /{dg} = T X /(dq) , /[dQ]' = T £ /(dq) 
where 



(dq) 



(2vr 



,£»-!■ 



Note that in Sect. 3 we keep the fermions as four-component objects, Tr(7 M 7^) = Ag^ 
eventhough the rest of the lagrangian is dimensionally continued. 

To simplify the reading of the forthcoming sections, we recall here the tensorial 
structure of the photon propagator D^ u (D~* = Dq ^ u + TL^) in the general case, and also 
in the hard thermal loop approximation. We follow closely the conventions of Ref. 
and write 

W»{P) = V^U L ( Po ,p) + V^U T ( Po ,p), (2.1) 

where P M = (po,p), p = |p|, p l = p l /p, and the subscripts L and T refer to longitudinal 
and transverse directions with respect to the vector p: 

= 7?« = o V% = S ij - p'pP 

n v + = - pi- - sr. (2.2) 

In the static limit (po = 0) the only non trivial components of D^ v are 



in the covariant gauge with gauge fixing parameter a. Note that ITl(0, p) = —Uqq(0,p) and 
n-r(0,p) = (l/2)Iljj(0,p). Please note that for simplicity most of our explicit calculations 
in Sects. 3 and 4 are in Feynman's gauge a = 1, but the gauge-independence of our main 
results (the propagator poles) will be demonstrated. 

The one-loop polarization tensor for ultrarelativistic gauge plasmas has been com- 
puted in Refs. [7[ |J. The dominant contribution in the high-temperature limit (i.e. for 
external momenta which are small compared to the temperature; e.g. po and p of or- 
der eT) is the hard thermal loop, which has the same structure for both abelian and 
non-abelian plasmasjl], £|. ||: 

nj2(P)^{- 9 ^„ + Po /f^}, (2.4) 
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where m 2 = e 2 T 2 /3 (with e 2 — > g 2 N for a SU(iV) gauge plasma), = (1, v), |v| = 1 and 
the angular integral / dQ runs over all the directions of v. We have indicated the order 
of perturbation theory as a superscript on the polarisation tensor, a convention that we 
shall systematically use througout. The structure ( |2.4| ) has a classical origin, as shown 
by thesimple kinetic derivation which is briefly reviewed in Appendix A. For soft P, 
Dq 1 ~ P 2 ~ e 2 T 2 ~ and the hard thermal loop must be included in the photon 
propagator. We denote this by *D )JLV : *D~l = Dq^ v + U*$. In the static limit, 



ng(o, P ) = -m 2 nl 2) (0,p) = 0, (2.5) 



so that 



* £>oo(0, p) = - 9 1 9 *D ij (0, p) = % . (2.6) 



9 , 9 *y V ' r/ 9 

pZ _|_ m z pZ 



2.2 Exact relations 



To derive eq. (|1 . 1|) , consider the plasma in the presence of weak static external sources 
with charge density p ext (x). The free energy in the presence of the sources is F = 
— lnTrexpj— j3(H + Hi)}, where Hi is the Hamiltonian describing the interaction 
between the gauge fields and the external sources: H i = J d 3 x p ext (x) A)(x). We assume 
that the average values of the gauge fields vanish in equilibrium, i.e. in the absence of 
sources. Then, to second order in p ext , the modification in the free energy reads 

F = F - l|(d P )p^( P ) A>o(0, P)p ext (~p) (2.7) 
1 f p ext (p) p ext (-p) 

= Fo+ 2i (dp) P 2 + n i( o, P ) ' 



where F is the free energy in the absence of sources, and D 00 (0, p) is the exact electrostatic 
propagator in equilibrium. By choosing p ext (x.) = q\5(x — xi) + ^^(x — X2), one extracts 
from the above equation the interaction energy of two isolated charges in the medium 
(r = Xi - x 2 ): 



V(r) = qiq 2 / (dp) — , (2.8) 

J p 2 + Tl L (0,p) 



which is eq. ( [1.11 ). 



One way to understand eq. (|2.7| ) is to recall that the free energy F[p ext ] is the gen- 
erating functional of connected Green's functions. Alternatively, one can express the free 
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energy in terms of the average gauge field Aq. This is achieved by performing the Legen- 
dre tansform -F'[A ] = F[p ext ] — J d 3 x p ext A . In F', the term quadratic in A involves 
the inverse propagator: 

F' = F + i|(dp)A (p) J D 00 1 (0,p)A (-p) + ... (2.9) 

= -If (dp) A (p) (p 2 + n L (0, p)) Ao(-p) + ... . 



Thus 



, , 1 5 2 F' 

n L (o,o) = 



(2.10) 

A =0 



V 5Al(p = 0) 

where V is the volume of the plasma. We note now that the chemical potential p enters 
the calculation of the partition function the same w= ay as Aq does, that is, it amounts 
to adding to the Hamiltonian a term — fi J d 3 :rp e (x) = — pp e (p = 0), where p e (x) is 
the charge density operator (in QED, p e (x) = V'^WV'W)- Thus, a change of A (p = 
0) is equivalent to a change of — p/e. Since F' and the pressure are related by the 
thermodynamic relation F 1 = —PV, eq. ( |2.10|) is the same as eq. ( |1.3|) . 

In order to establish eq. ( |1.4j ), we rely on the exact Dyson-Schwinger equations and 
on the Ward identities. This will also provide us with an alternative proof of eq. ( |1.3| ). 
Notice first that gauge symmetry ensures the transversality of the polarisation tensor, 

P^(p ,p) = 0, (2.11) 

so that 

pTT„(0,p) = 0, (2.12) 

which implies ITo(0,p) = and the transversality of ITj(0,p) in the spatial indices. 
For spinor QED, the relevant Dyson-Schwinger equation reads 

IV(po,p) = -e 2 f {dK} Tr( 7 ^(P + K)Y V {P + K, K)S(K)) , (2.13) 

where S is the full fermion propagator, T the full vertex, and the trace is over spinor 
indices. In the limit po = 0, p — > 0, eq. ( |2.13| ) combined with the Ward identity T U (K, K) = 
dS-\K)/dK v gives 

n M ,(po = 0,p ^ 0) = -e 2 Tr 7m j {dK} . (2.14) 

In the imaginary-time formalism, the variable k° = iuj n + p e takes discrete values only. 
Then, the derivative with respect to k° in the above equation is meant to be done on the 
analytic continuation of S(k ). 
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In the calculation of II(oo) from the r.h.s. of eq. ( |2.14|) , one can replace d/dk 
d/dfi e ; then, by noticing that 



1i lp J{dK}S{K) = 3p = g pQ p e (2.15) 



is the charge density in equilibrium, and that p e = dP/dp e , we recover the identity ( |1.3|) . 
As for the components, they vanish after an integration by parts, in agreement with 
eq. (0)- 

A similar discussion may be carried on for SQED. The corresponding Dyson- Schwinger 
equation is illustrated in Fig. 1. To save writing, we consider directly the limit po = 0,p —>■ 
in these diagrams, and denote U pu = U pu (po = 0,p — > 0) in the remaining of this section. 
It is convenient to combine the two diagrams in Figs, l.a and l.b by writing 

K» + = -2g, v e 2 J [dK] S(K) - 2e 2 J [dK] K p T V {K, K) S 2 (K) (2.16) 

where S(K) = l/(— K 2 + H(K)) is the exact scalar propagator, and T V (K, Q) is the vertex 
with one photon and two scalar external lines, with K (Q) denoting the momentum carried 
by theincoming (outgoing) charged particle. We use the same Ward identity as in spinor 
QED to rewrite eq. ( |2.16|) as 

+ IT^ = -2e 2 J [dK] —(K,S(K)) . (2.17) 



After an integration by parts (with the assumption that kiS(K) — > as |k| — ^ oo), we 
[ a . + n b . 



obtain II"- + H|- = 0. As for the electric piece, this is rewritten by introducing a small 



chemical potential for the charged particles, so that k° = iu n + /i e and 

ng + 14 = -2e 2 — J [dK] (k S(K)) . (2.18) 

The contribution of the remaining diagrams, Figs, l.c and l.d, is evaluated as 

+ 11%, = 2e 2 J [dK] I [dQ] {g° p S{K) D pa {K - Q) S(Q) T pv (K - Q, 0, K, Q) (2.19) 

- 2S 2 (K) T U (K, K) S(Q) D m {K - Q) T CT (Q, K) 
where V pv is the vertex between two photons and two charged scalars. This vertex satisfies 

V pv {K - Q, 0, K, Q) = -e UL + JL^j r p (K, Q) . (2.20) 

By using the Ward identities above, one obtains II? • + = and 

W m + 14 = -2e 3 — J [dK] J [dQ] (S(K) S(Q) D 0p (K - Q) T P (K, Q)) . (2.21) 
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We recognize in the r.h.s. of eqs. ( |2.18| ) and ( [2.2 1| ) the average electric charge density 



p e = i(4>^d 4> — (9 O 0^)0) — 2eA (p ji (f) expressed in terms of exact propagators and vertices. 
Thus 

rioo = n° + n* + n^ + n^ = - e 2 |^ = -e 2 1| . (2.22) 

It is interesting to observe that the vanishing of Uij is obtained by independent cancel- 
lations of the diagrams l.a and l.b among themselves, and of the diagrams l.c and l.d 
among themselves. 

Note that the manipulations above are formal in the sense that we have left aside 
the question of the regularization of UV divergences, as well as possible IR problems. We 
will not address these questions in general here but in the following sections we shall go 
through an explicit perturbative verification. 



3 Perturbative QED 
3.1 The static- infrared limit 

We now begin our explicit calculations by considering in this subsection the object 
n M ;,(0,p — > 0) up to fifth order (e 5 ). The explicit computations here have been performed 
in the Feynman gauge. At one loop, a standard calculation gives 

Tr (7^7*40 



ng)(o,j>-o) = wj{dK) 



-2 



-4W)*J{dK} 9 '» K J*"*" , (3.1) 



with If[ = j p K p '. As usual, /i is the mass scale introduced by dimensional regularisa- 
tion, and the gauge coupling e is dimensionless. Within dimensional regularisation, an 
integration by parts in the k integral gives 

{dK] Ki = — 2 — y-ffT- (3 - 2) 



Hence one deduces from (|3.1| ), 

ng } (0,p^0)=0, (3.3) 

and 

,2_ TT(2)/ n „ _ A/ „..^2 fe) n\fi dK } 



m 



nS(0,p^0) = -4(e//) 2 (2 -D)J (3.4) 

2 T 2 

= (3-5) 
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where the limit D — > 4 has been taken in the second line. In the rest of this paper, this 
limit will often be taken implicitly in our final results. The fermionic integral occurring 
above is a special case of 

fn = / = (2 2n+1 ^ " 1) (3-6) 



(Q 



where 



[dQ] 

(Q 2 ) n 

= T D ~ 2n 2 \ % C(2n + 1 — D) T P + , (3.7) 

(2tt) 2 " r(n) V 7 V 2 / ' v ; 

and ((x) is Riemann's ^eta-function. The quickest way to evaluate the b n is to first 
integrate over the momenta and then perform the frequency sum. The integral f n is then 
obtained using the following trick: consider the sum b n + f n and rescale the momenta 
there by a factor of 2, giving b n + f n = 2 2n+1 ~ D b n , from which the quoted result follows. 
Note that, because of the dimensional continuation, the Matsubara mode go — does not 
contribute to the integral in eq. (|3~?1 ). 

The next contribution is of order e 4 and comes from the two loop diagrams shown 
in Fig. 2. Each of these diagrams is quite complicated, even in the static IR limit, but 
their sum is remarkably simple and is given by 

n<2(o )P - o) = A{e^)\D - 2)(rr> ) + n^) (3.8) 

with 

njf;) = J{dKdR} K » R ^ R » , (3.9) 

= (61 - h) J {dK} (|^ - . (3.10) 

Note that eq. ( |3.9|) vanishes at zero chemical potential (at nonzero chemical potential it 
only contributes to il o)- Using integration by parts, as in (|3.2| ), one obtains 

4 4) (0,p^0) = 0, (3.11) 

and 

Ilgtop-O) = iM 4 (D-2)(D-4)(b 1 -f 1 )f 2 = (3.12) 

Notice the cancellation between the UV divergencies of eqs. ( |3.9| ) and ( p.!0| ) in their sum 
( |3.12|) (the quantity (D — 4)/ 2 being finite for D — > 4). This cancellation is necessary 
since the contribution from the UV counterterm diagrams mutually cancel (see later). 
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Consider now the behaviour of the magnetostatic polarization tensor for small, 
but non- vanishing, momenta. We shall argue here that Uij(0,p — > 0) = 0(p 2 ) to all 
orders in perturbation theory. To see this, consider an arbitrary graph contributing to 
Uij(0,p). Label the photon lines in the graph by a set of independent loop momenta so 
that the external momentum p only flows along the fermion lines. Since in imaginary 
time the fermionic propagators are infrared safe, they can be expanded with respect to 
the soft momentum p. Then the external momentum appears only in the numerator and 
rotational symmetry ensures that the terms with odd powers of p vanish. Hence, IIy(0, p) 
is analytic in p 2 for small p, so that n^(0, p — > 0) = 0(1) + 0(p 2 ). By also using eq. (|1.4|) , 
we conclude that Uij(0,p) vanishes at least as p 2 as p — > 0. 

The proof in the last paragraph is rigorous only when applied to diagrams without 
photon self-energy insertions since these latter cause power-like infrared divergences along 
internal lines due to the nonvanishing of IToo(0,0), thus making the arguments formal. 
For diagrams with self-energy insertions, the arguments of the last paragraph must be 
applied not to individual diagrams but to sums of similar diagrams which result in a new 
effective graph with dressed photon propagators. The conclusion then is as before, namely 
Uij(0,p — > 0) = 0(p 2 ). Finally, note that the proof is independent of the fermion mass. 

Beyond fourth order we begin to get contributions nonanalytic in e 2 . The e 5 term 
comes from dressing the zero mode of the photon propagators, as shown in Figs. 3 (cf. 
fl5H). The nonzero modes of the photon line (like the modes of the fermion lines) are cut 
off in the infrared by the scale T and dressing those modes just gives the usual perturbative 
corrections (which are analytic with respect to e 2 ). Consider now a static internal photon 
line: inserting the electric polarisation tensor along this line causes infrared divergences 
which can be summed up, with the result that the bare electrostatic propagator gets 
replaced by *D 00 (0, q) = — l/(q 2 +m 2 ), m 2 = e 2 T 2 /3 (recall eq. (p76|) ) . On the other hand, 
insertions of the static magnetic polarisation tensor give only perturbative corrections, as 
one can verify from power counting by using n« (0, q) = 0{q 2 ) (also recall that n i(0, q) = 
to all orders so that there are no corresponding insertions). 

The fifth order contribution is then obtained from the diagrams in Fig. 3, in which 
the internal photon line is static and is associated with the dressed propagator *D^ U (0, q), 



eq. ( |2.6| ). It is then convenient to write 



>iV(o,q) = AUo,q)+ f 9 2 T\ 

q (q + m z ) 



= D ,^(0,q)+*d^(0,q), (3.13) 
and to observe that it is the piece *d which is responsible for the e 5 contribution. In 
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fact, to order e 5 , a further approximation can be done. This consists in neglecting the q 
dependence along the fermion lines. Then, the q integral decouples and is easily evaluated: 
it is of order / d 3 qm 2 /q 2 (q 2 + m 2 ) ~ m ~ e, and, since one has a factor of e 4 from the 
four vertices, the net result is of order e 5 . The forgotten g-dependence along the fermion 
lines only gives a sub leading contribution because the fermion lines are IR safe, so that 
one may expand out the q dependence and do the usual power counting. 

Thus the exact fifth order contribution comes from the sum of the diagrams in Fig. 
3 when *d(Q) is used for the photon lines and the Q dependence along the fermion lines 
is ignored. We get: 

n{5(o,o) = ^ m "r / J (3-14) 

S, u = (4K,K u -2g, u K 2 )(K 2 + 2k 2 ) + K\2g, g u0 -g^) 

-%k K 2 {g^K v + g u0 Kj + lGK^K 2 + k 2 ) . (3.15) 

Again using manipulations as in ( |3 .2D gives 

nS?(0,P-0) = 0, (3.16) 



where 



and 



ngCO.p-O) = A(e^m D - 3 T(D-2)(D-A)f 2 J 



(dx) 



x 2 (x 2 + 1) 



-e 5 T 2 



4tt 3 V3 ' 



(3.17) 



In summary we have shown for massless QED in the Feynman gauge at temperature T 
and zero chemical potential /i e , 

Ua{0,p -> 0) = 0(p 2 ) + 0(e 6 T 2 ) . (3.18) 

and 

n 00 (0. P 0) = g - ^ + ^) + 0{ eV) . (3,9) 



A scan of the computations shows that Eq.( |3.18 ) has also been verified for nonzero chem- 



ical potential. The result ( |3.19| ) — which is the left-hand-side of the identity eq. flO] 



is in accordance with the r.h.s. of as given in p|. Also, since the right-hand-side of ([O 
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is gauge- independent, this implies the gauge- independence of ( |3.19| ) though we have done 
our calculations in the simpler Feynman gauge. 

The result eq. ( |3.19| ) has turned out to be UV finite eventhough we seem to have 
ignored the UV renormalisation. However recall first that zero-temperature counterterms 
suffice to render the theory UV finite, and that gauge invariance ensures that the vacuum 
(i.e. zero-temperature) polarization tensor vanishes in the zero-momentum limit. Then, 
in our calculation, no counterterms were needed for the subdiagrams of the two loop 
polarization tensor for the following reasons: i) The order e 4 vertex and fermion wave- 
function counterterm diagrams cancel against each other because of the Ward identity 
Z\ = Z 2 ; ii) since we are working in the massless limit, there is no mass counterterm; 
iii) there is no order e 5 counterterm diagram. On the other hand, UV divergences will 
make their appearance when the polarization tensor will be considered for non-vanishing 
momenta in the next section. Similar arguments apply to the SQED calculation in Sect. 4, 
and so we shall not repeat them there. 

3.2 Screening masses 

Having computed n M „(0, p — > 0), we are ready to determine the screening masses. To one- 
loop order, the electrostatic propagator *Doo(0, p) = —1/ (p 2 +m 2 ) has a simple pole in the 
upper half of the complex p-plane, occuring at p = im. Accordingly, m 2 D = m 2 + 0(e A T 2 ) 
(recall eq. (|1.2|) ). The higher-order corrections that we have calculated modify the position 
of this pole, without changing the analytic structure of D O o(0, p). Accordingly, we may 
define the electric screening mass as the solution of the equation 

p 2 -U 00 (0,p) = (3.20) 

for p ~ eT. Since we will only do the calculations up to order e 5 , we need the expansion 
of the one and two loop static polarisation tensors up to the following orders, 

nS(0,p~eT) = -e¥(a 20+ ^ + ^ + ^) + O(eV/^) (3.21) 
n$(0,P~eT) = e¥(a 4O + ^) + 0(eV), (3.22) 

together with nj§(0,0). 

From the last subsection we have a 2 o = 1/3 and a 40 = l/8n 2 . Furthermore a 2 i = 
&23 = Q-41 = for the reasons explained after eq. ( |3.12| ). We now calculate a 2 2 starting 
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from 



Ilg^O.p) = 2(e//) 2 J \dK}5 Po>0 ^ + 



Ak 2 - p 2 
K 2 (K + P) 2 



(3.23) 



Since ko is discrete and of order T, one may expand 



+ P) 2 K 2 



2k-p + p 2 (2k-p + p 2 ) 2 
1 H r ,„ + - 



K 2 



if 4 



(3.24) 



Then using (|3.24 ) in (|3.23| ) and simplifying , we get 



ng(o,p) = n$(o,o) - (D - 2)/ 2 + 0(eV/T 2 ) . 



(3.25) 



The order e 2 p 2 term in (|3.25|) is UV divergent as D — > 4 and this divergence is cancelled 
by the photon wave-function renormalisation counterterm 5Z 3 p 2 , with 5Z 3 = — e 2 /127r 2 e. 
Then the renormalized value of reads 



a R 



4/i 2e (l 



-h 



1 



7 - 1 + In 



Afi 2 

J^2 



7T 



3 '* 127r 2 e 12tt 2 
We summarize the above results by writing (p ~ eT) 

n 00 (0,p) = -m 2 - af 2 eV + 11$ (0, 0) + ng(0, 0) + 0(e 6 T 2 ) 

The solution p 2 = —m 2 D of eq. ( |3.20| ) to fifth order is therefore 



(3.26) 



(3.27) 



m 



D 



T — 



8tt 2 36tt 2 



7 - 1 + In 



V 
ttT 2 



+ 



+ 0(e 6 T 2 " 



(3.28) 



Here, e = e(/z) is the running coupling constant defined by the minimal subtraction 
scheme. Since this satisfies de/dln/x = e 3 /127r 2 , it is clear that the r.h.s. of eq. ( p.28|) is 
independent of the renormalization scale /i, to order e 5 . 

Eq. ( p.28| ) is our result for the electric screening mass of massless QED at tempera- 
ture T and zero chemical potential. Up to order e 4 this coincides with the result of Ref. ||. 
As for the magnetic screening mass, this vanishes as expected, since IIy(0, p — > 0) = 0(p 2 ). 

Let us now establish the gauge-independence of our result ( p.28|) . The constants am 
come from the one- loop diagram and so are manifestly gauge- independent. The constant 
040 is gauge-independent because of the relation (|1.3Q while the vanishing of 041 is a gauge- 
independent statement since the arguments following eq. ( |3.12j ) make no reference to any 
gauge-choice. Hence the electric screening mass (and similarly the vanishing magnetic 
screening mass) to fifth order as given by eq. (|3.28| ) is gauge-independent. 
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4 Perturbative Scalar QED 



In this section we consider the electromagnetic interactions of a charged scalar field </> 
described by the Lagrangian (D^ = <9 M + ieA^) 

£ = - ^(3 ■ Af + (D^(D^) - -^<pf. (4.1) 

Since we are interested only in the effects of the electromagnetic interactions, we shall 
ignore the self coupling of the complex field 0, i.e we assume A — > in what follows. 

In this section, we shall compute the correction of order e 4 T 2 to the Debye mass, 
and we shall verify that 11^(0, p) = 0(p 2 ) as p — > 0, to the same order. In order to 
avoid double counting in higher order calculations, we shall refer to the skeleton diagrams 
displayed in Fig. 1. The corrections to U^ u at various orders will be obtained by expanding 
to the appropriate order the exact propagators or vertices in these diagrams. As explained 
before, in this procedure one must keep the thermal masses on the static propagators, 
while non-static propagators can be perturbatively expanded. 



4.1 Leading order results: hard thermal loops 

In leading order, the self-energies for the photon and the scalar particle are obtained from 
the 1-loop diagrams in Figs. 4 and 5 respectively. For soft external momenta (that is, 
for po and p of order eT) , the dominant contributions to these diagrams come from loop 
momenta of the order of T: these are the hard thermal loops]!], ||. The photon hard 
thermal loop was already presented in Sect. 2 (eqs. (|2.4j )- (|2T^ )). The hard thermal loop 
for the scalar self-energy reduces to a (gauge-independent) local mass term: 

S ( 2 )(P) = {D-1) (e^) 2 /[dQ] So{Q) = — ^— = M 2 , (4.2) 

where S (Q) = -1/Q 2 . We define the propagator *S by *S~ 1 = Sq 1 + £( 2 ). In the static 
limit, 

Generally, the one-loop results cannot be trusted beyond the hard thermal loop 
approximation. This is so because, beyond leading order, soft loop momenta start to con- 
tribute and the corresponding propagators must include the hard thermal loops. Consider 
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for example the small momentum behavior of the static one-loop polarisation tensor. It 
is shown in Appendix B that 

Ill L (0,p) = ±e 2 pT + O(ey) (4.4) 
(see eq. ( |B.2| )). A similar behaviour is obtained for the scalar one- loop self-energy: 

£n,(0, p) = M 2 — \?pT + 0(e 2 p 2 ). (4.5) 

These non-analytic (in p 2 ) contributions arise from the static internal modes. We shall 
verify in the next subsection that these terms disappear once thermal masses are included 
in static internal lines. In QED, such a problem does not occur because the corresponding 
one-loop diagram for IL^ has only fermionic internal lines. Note that for the component 
IIoo no resummation is needed to get the leading low momentum behavior to one-loop 
order (see eq. (|B.6| )): because of the vector structure of the electromagnetic interaction, 
the mode n = does not contribute to the Matsubara sum in the contribution of the 
diagram 4.b to U 00 (0,p); as for the diagram 4. a, it is momentum independent. 



4.2 Next-to-leading order: ring summation 

In this subsection, we consider the consequences of the resummation on the one-loop 
diagrams of Figs. 4 and 5. We isolate the static Matsubara mode, since this is the 
only one which is concerned by the resummation, and we replace the bare propagators 
by the propagators *D^ V and *S obtained in the hard thermal loop approximation. We 
thus obtain the dressed one-loop diagrams of Figs. 6 and 7, whose contributions will be 
denoted by * 11^(0, p) and *S(0,p), respectively. 

The tadpole diagram in Fig. 6. a gives 

= -2g, v (e//) 2 T / (dk) (_L_ - i.) . (4.6) 



The second term inside the brakets in eq.( |4.6| ) substracts that contribution of the static 
mode which has already been included in the hard thermal loop, eq. (|2.4j) . Note that this 
term vanishes in dimensional regularisation so we could have as well omitted it. However it 
is better to keep it here. This will allow us to show that the final results, when appropriate 
contributions of a given order are added, are both ultraviolet and infrared finite, even in 
the absence of regularisation. Letting D — > 4 in eq. (fOl), we get 



*™ = 9^ 2 MT = 9^ 3 T 2 (4J) 

^ 2vr 4vr y j 
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The diagram in Fig. 6.b gives a non trivial contribution, of order e 3 , only to the magnetic 
piece lljj(0,p). As p — > 0, 

U*(0, 0) = -4 (e^) 2 T / (dk) k 2 - ^) (4.8) 

= A e 2 MT = — e 3 T 2 , 
27T 4tt 



By adding together the contributions Q4.7|) and (|4.8j ), we obtain the total contribution of 
order e 3 to the zero-momentum limit of the polarisation tensor 0: 

p 2 MT p 3 T 2 

nS } (o,o) = i^i = nf (o, P - o) = o. (4.9) 

The above correction to n o(0, 0) can be understood as a classical correction, as discussed 
at the end of Appendix A. 

Consider now the momentum dependence of *Il Atl/ (0,p). The electric piece is inde- 
pendent of p, since it is entirely given by the tadpole of Fig. 6. a (recall eq. ( [4.7| )). For the 
magnetic piece, both diagrams in Fig. 6 contribute, and give 

*n«(o,PH WT/(dk) {^- (k , + X (k + + P) p 2 )2 + M2) }. (4-10) 

The presence of the mass M ~ eT in the denominators allows an expansion of the last 
one with respect to p. After integration over k, only the terms even in p survive. The 
small momentum expansion of Q4.10Q is therefore 

•n«(0, p) = ^ ^ {l + c a (p/M) 2 + c 2 (p/M)* + ...}, (4.11) 

where the c^'s are constant coefficients. Note that there is no term linear in p, contrary 
to the pure one-loop result of eq. ( B.l| ). As p — > 0, the leading term is proportional to 
e 2 (T/M)p 2 ~ ep 2 . For p ~ eT, all the terms in the r.h.s. are of order e 3 T 2 , and the 
integral in ( f4.10|) should be computed exactly, with the following result: 

e 2 MT f 4M 2 + v 2 V 1 

*n«(0, p) = P arctan — - 1 ■ (4.12) 

v 2vr \ 2pM 2M J K J 

For small momenta p ^ eT, this equation gives the leading infrared behaviour of the static 
magnetic polarization operator (the non-static modes in the one-loop diagrams in Fig. 4 
contribute only to order e 2 p 2 ). For large momenta, p/M 3> 1, *IIjj(0, p) — > e 2 pT/8, as for 
the undressed one-loop result of eq. Q4.4|) . 
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A similar discussion applies to *E(0, p), the correction to the static scalar self- 
energy given by the dressed one-loop diagrams of Fig. 7 (which, we recall, involve only 
the internal mode with zero frequency). A straightforward computation gives 

r f 1 1 D — ^ 2 n 2 — M 2 1 

• S (0,p) = (e,fT/(d q) |^ + + - ^ ( J p)2 + M2 }. (4.13) 

in the Feynman gauge. The g-integral can be performed readily, and for (D — 1) — > 3, 
one gets 

» (2(M 2 -P\ rctanjL _ !H±i£\ , (4 . 14) 

V P; 4tt \ Mp M My y 1 

which admits the following small-momentum expansion: 

•2(0, p) = t^L _ I^Mf + 0(P VM*)} . (4.15) 



For high momenta, p/M ^> 1, we recover the linear behaviour in p as in eq. ( |4.5| ). For 
p ^ eT, eq. (|4.14| ) gives the dominant non-trivial momentum behaviour of the scalar 
self-energy. 



4.3 Order e 4 : diagrams l.a and l.b 

Contributions of order e 4 T 2 arise from two-loop diagrams in which the static propagators 
are dressed by thermal masses. Because of the resummation involved in the thermal 
masses, parts of these two-loop diagrams have already been included in the order e 3 
calculation. In order to avoid double counting, we refer to the skeleton diagrams of Fig. 
1. We shall give details only for the tadpole diagram, Fig. l.a: 

U% = -2g, v (e//) 2 J [dK] S(K) = g, u Yi a (4.16) 

where S(K) = l/(—K 2 + S(i^)) is the exact scalar propagator, with S(i^) the exact 
self-energy. In these expressions K = (2iirnT, k). When n ^ 0, S represents a small 
perturbative correction and only the second term in the expansion 

SW = -^(l + ^ + -) (4.17) 

is in fact needed to evaluate ( |4.16| ) up to order e 4 . When n = 0, infrared divergences 
render the expansion above meaningless. It is then convenient to expand about the 
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massive propagator *S(0, k) (eq. (|4.3|)): 



S(0, k) 



k 2 + M 2 + (S(0, k) — M 2 ) 



(4.18) 



+ 



1 



1 



E(0, k) — M 2 



+ ■ 



k 2 Vk 2 + M 2 kV (k 2 + M 2 ) 2 
When the expansions ( |4.17| ) and ( |4.18|) are used in eq. (|4.16| ), the following result is ob- 
tained 



rr 



-2 (e//) 2 J[dK]S (K) 

-2 (e//) 2 T J (dk) [*S(0, k) - S o {0, k)] 

+2(e^) 2 /[dATS 2 (K)S(K) 

+2 (e//) 2 T / (dk) *S 2 (0, k) [E(0, k) - M 



(4.19) 



where the neglected terms are, at least, of order e 5 T 2 . The self-energy entering the r.h.s. is 
the one-loop self-energy, that is, it is obtained from the diagrams in Fig. 5 or, if necessary, 
from the dressed diagrams of Fig. 7 (see below). 

Consider now the different terms in the right hand side of eq. ( 4.19 ). We have 
already evaluated the first two integrals giving respectively the contributions of order e 2 
and e 3 . For the third integral, it is enough (since k° ~ T) to use the one-loop expression 
of the scalar self-energy S, i.e. 

(4.20) 



(4.21) 



2(e^) 2 J[dK]'S 2 (K)i: 1L (K), 
where (see Fig. 5) 

E 1L (K) = M 2 + 2 (e//) 2 K 2 J [dQ] S (Q) S (K + Q), 
in the Feynman gauge. 



The evaluation of the last term in eq. (|4.19|) is more involved. It is again necessary 
to separate the static (u m = 0) and non-static {uj m ^ 0) modes in the one-loop diagrams 
giving S(0,k) (see Figs. 5 and 7); here, u m denotes the Matsubara frequency inside the 
loop. For the non static modes, bare propagators can be used, as in Fig. 5, and we 
recover the m ^ piece of the one- loop self-energy from eq. ( |4.21| ). The corresponding 
contribution to IP reads 



-4(e//) 4 T|(dk) 



2(e//) z T /(dk)*S 2 (0,k) 
k 2 



(k 2 + M 2 y 



[dQ]' 



Si L (0,k) 
1 



1 



(4.22) 



+ q 2 + (q + k) J 
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It is not difficult to see that the above contribution is of order e 4 MT ~ e 5 T 2 (in D = 4) 
and therefore can be ignored in our present computation of the order e 4 . In fact, eq. ( 4.22] ) 
is precisely of the type already encountered in Sect. 3.1, in computing the contribution 
of order e 5 to 11^(0,0) for spinor QED. In particular, its ultraviolet singularity in the 
limit D — > 4 is harmless, since it will be compensated by similar contributions arising 
from other mixed two-loop diagrams (as happens, e.g., in eq. (|3.17|) ). (We characterize as 
"mixed" any two-loop graph where one of the internal frequencies is non-vanishing, while 
the other one is zero.) Mixed graphs do not contribute to 11^(0, p) to order e 4 , and will 
be neglected in what follows. (This applies, in particular, to the mixed graph included in 
eq. CT .) 

Consider now the remaining contribution to eq. (|4.19|) , that is, 



2(en e ) T /(dk)*5 2 (0,k) *£(0,k)-ilT 



m=0 



(4.23) 



This corresponds to two-loop diagrams where both the internal frequencies are zero, so 
that the corresponding propagators are dressed by the hard thermal loops. The integral 
( f4.23|) is explicitly written in eq. (|4 . 3 2|) below. 

We turn now to the diagram l.b, which gives 

n^(0) = -2 (e//) J [dK] k, Y V {K, K) S 2 (K) . (4.24) 

To one-loop order, and for Q = K, one readily gets T U (K, K) = 2e[fK v + 5Tl L (K,K), 
with 

8Tf{K, K) = -4 (e//) 3 j [dQ] S (Q)S (K + Q) {k, + (K, + QJ&S^K + Q)} (4.25) 

in Feynman gauge. For k = iu n ^ 0, the one-loop expressions for the vertex and 
the scalar self-energy, as given by (|4.21| ) and ( [4.25D , are sufficient to get the order e 4 
contribution to eq. ( |4.24j ): 



nit 6) (o) 



-2 (e^) / [dK]' k, 5Tl L {K, K) S 2 (K) 



(4.26) 



+8 (e//) 2 J [dK]' kyk v E 1L (K) S 3 (K). 



In the calculation of n o we only have to consider non vanishing Matsubara frequencies 
ko 7^ 0, so that Hocf^ is completely determined by eq. ( (4.26| ) above. The static mode only 
contributes to the magnetic sector, and there the resummation of the thermal masses is 
again necessary. In this case, we expand as follows (compare eq. (|4.19| )): 



n=0 



2 (e/i e ) T / (dk) ki 5Ti(k, k) *^ 2 (0, k) 



(4.27) 



+8(e/i e )'T /(dk)k 2 *^(0,k) S(0,k)-M 
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Let u m denote the internal Matsubara frequency in the one-loop diagrams contributing 
to <5rj(k, k) or to X(0,k). The non-static modes (u m ^ 0) give rise to mixed two-loop 
graphs which contribute to eq. ( |4.27| ) only to order e 5 ; they will be ignored in the present 
calculation. For the static mode (u> m = 0), the resummation of the thermal masses is 
compulsory, so that the corresponding self-energy *£(0, k) and vertex-correction *5Ti(k, k) 
are determined from dressed one-loop diagrams. The corresponding result for the self- 
energy has been given in eq. ( |4.13|) . For the vertex function, one obtains similarly 



[fci**ri(k,k) 



\m=0 



= -4 (e//) 3 T J (dq) S o (0, q) *S(0, k + q) 

{k 2 - k • (k + q)(k 2 - M 2 ) *S(0, k + q)} 

<9*£(0,k) 



(4.28) 



(e//) h 



m=0 



As shown by the last line above, this result is consistent with the Ward identity and 
with the expression (^T3|) of *S(0, k). After combining eqs. fl4~27l) , (^T3|) and ( ^28|) , one 
obtains the result displayed in eq. ( [4.39|) below. 



We now summarize the results obtained in this section. We have 



where 



n( 4a ) = A 1 + A 2 + A 3 



Ai = 2 (e/i e ) 2 M 2 J [dK]' S$(K) 



(4.29) 



(4.30) 



A 2 = -4 (e/i e ) 4 J [dK]' J [dQ]' S (K)S (K + Q)S (Q) 



(4.31) 



A 3 = 2(e^) 4 T 2 |(dk) |(dq) 



1 



(k 2 + M 2 ) 2 
2 2 k 2 -M 2 1 
q 2 + m 2 + q 2 + M 2 _ q 2 ~ q 2 (q + k) 2 + M 2 J ' 



1 



1 



(4.32) 



The first two terms, A\ and A 2 , are obtained from eqs. (|4.20| ) and ( 4.21 ). In writing them, 
we have omitted mixed terms with uj n ^ and u m = 0, since they contribute only to 
higher orders. The term A 3 follows immediately from eqs. ( |1.23| ) and ([4.13|) . 

Similarly 



2 

00 > 



(4.33) 
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with 



i4 = 8(e//) 2 M 2 |[dK] k 2 S 3 (k), 

B 2 = -12 (e//) 4 J [dK] J [dQ] k 2 Q S 2 (K) S (K + Q) S (Q). 
Finally, for the magnetic contribution we set 



(4.34) 



(4.35) 



(4.36) 



with 



Bl = -2{e^) 2 M 2 J[dK]k t ^-Sl{K) 



(4.37) 



B 2 = 4 (e//) 4 J [dK]' J [dQ]' S (Q) k t — [S (K)S (K + Q)] , (4.38) 



B\ = -2 (e^) 4 T 2 /(dk) j(dq)k 







1 



+ 



dtt { (k 2 + M 2 ) 2 
2 k 2 — M 2 



q2 + m 2 q 2 + M 2 q 2 q 2 (q + k) 2 + M 2 



(4.39) 



One can verify that the magnetic contributions from diagrams (a) and (b) compen- 
sate, as expected. To see this, note that we can write 



nf) = (l>-i) /(dk)/(k) 



and 



nl 4fe) (o) = /(dk)^A /(k) . 



(4.40) 



(4.41) 



The sum of the two expressions above vanishes after an integration by parts (allowed by 
dimensional regularisation). One can also verify that u\f a \o, p) + 1X^(0, p) vanishes as 
p 2 for p —>■ 0. To do this, it is sufficient to study the contribution of the static modes for 
Po = but non- vanishing p (that is, to replace Bf t , eq. ( f4.39| ), with the corresponding 
contribution for p ^ 0). Then, it may be verified that, because of the resummation of 
thermal masses, the relevant expression admits a well-defined expansion in powers of p 2 . 
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4.4 Order e 4 : diagrams l.c and l.d 

With bare vertices, and zero external momenta, the contributions coming from the dia- 
grams l.c and l.d are, respectively, 

n^(O) = -4 (e//) 4 J [dK] J [dQ] D, U (Q)S(K)S(K + Q) , (4.42) 

and 

n^(0) = -8 (e//) 4 J [dK] J [dQ] D^(Q) (2K° + Q a ) K v S 2 (K)S(K + Q) , (4.43) 

with ko = iuj n and go = iw m . (The vertex corrections in these diagrams play no role up to 
the order e 4 .) If u n ^ and u rn ^ 0, bare propagators can be used in these expressions in 
order to obtain the e 4 contribution. If u n = uj m = 0, the propagators *S and *D should 
be used instead. Finally, the mixed graph where u n = and u m ^ 0, or vice-versa, do 
not contribute to order e 4 . In all cases of interest, the calculation is straightforward and 
leads to the following results: 

na c) (0) = C 1 + Co 2 o> (4-44) 
C l m = 4 (e//) 4 J [dK] (1 - S n0 ) J [dQ]' S (K) S (K + Q) S (Q), (4.45) 

C 2 = 4 (e^) 4 T 2 /(dk) /(dq) ^ ^ { ^ + AP ' ^ 

n!- c) (0) = Cl + Cl (4.47) 
C\ t = -4(D - 1) (e//) 4 J [dK]' J [dQ]' S (K) S (K + Q) S (Q), (4.48) 

d = -4P - 1) (e^) 4 /(dk) /(dq) £ ^ + ^ + ( 4 ' 49 ) 
nS d) (0) = Doo = 12 (e//) 4 / [dK] / [dQ] k\ S 2 (K) S (K + Q) S (Q). (4.50) 

nr } (0) = ^ + ^, (4.51) 
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D\ = 8 (e//) 4 J [dK]' J [dQ]' [k ■ (2k + q)] S 2 (K)S (K + Q)S (Q), (4.52) 



Dl = 8(e^ T 2 /(dk) |(dq) 



1 k • (2k + q) 1 

q 2 (k + q) 2 + M 2 k 2 + M 2 ' 



(4.53) 



It can be easily verified that the sum of the magnetic contributions of diagrams (c) 
and (d) is a total derivative with respect to k and therefore vanishes upon integration. 
For instance 



Cl + Dl = -4 (e//) 4 T 2 J (dk) J (dq) 1 



1 



A-' 



<9F 



(4.54) 



^k 2 + M 2 (k + q) 2 + M 2 / 



0. 



and similarly C\ + D\ = 0. One can also verify that the function llffi(0,p) + 11^(0, p) 
is analytic in p 2 as p — > 0. 

4.5 Electric mass to order e 4 



By adding together the results of sections 3.3 and 3.4, we are now able to obtain the cor- 
rection of order e 4 T 2 to n 00 (0, 0). The relevant equations are (|OD|)-(^2|, (g3J)-(g3J), 
(fT45])-([l]4(|) and ( |4.50| ), which express the relevant contributions of the four diagrams in 
Fig. 1. A simple inspection of these equations shows that the following compensations 



arise: 



M + Co 1 , 



no 



0. 



00 



0. 



(4.55) 



That is, the only non-trivial contribution of the non-static two-loops diagrams is that 
given by A\ + Bq , eqs. ( ^.30| ) and ( [4.34| ). It is proportional to the scalar thermal mass 
squared M 2 : 

A l + B l m = 2(e^) 2 M 2 J[dK]'S 2 (K) (l+4k 2 S (K)) (4.56) 
= 2M 2 (D-4)(e/i £ ) 2 6 2j 

where the second line follows from the first one after an integration by parts (see eq. ( |3.7| )). 
By evaluating b 2 from eq. (|3T7|), we obtain 



oo 



eMY 
~2W 



4:71 



(4.57) 
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for the contribution of the non-static modes. The contribution to ripo (0, 0) which arises 
after resumming thermal masses in purely static two-loops diagrams involves the terms 
A 3 and Cq (eqs. ( |4.32| ) and ( |4.46 )), whose sum is computed in Appendix B: 



^2 eT \ fm-M m + 2M\ , 4r . 

A 3 + CL = - h4 n (4.58) 

3 00 I 4tt / V M 2M J K J 



e 2 T 



2 




\ 47T J 

By adding eqs. ( [4.571 ) an d (|4.58|) , we derive finally 



1 + 4 In 1 




a), . f e 2 T\ 2 (m m + 2M\ 

nS»(0,0) = -(-) (- + 4.0^^) (4.59) 

- -©"(**«■('♦*))■ 

No UV divergencies have been encountered in the derivation of this result: all the terms 
which survive the compensations ( |4.55[ ) are UV finite. Note that the first line vanishes if 



m = 0, i.e. if we "forget" to resum the static longitudinal photon lines. 

The result ( f4.59|) represents the main ingredient for computing the correction of 



order e 4 to the Debye screening mass m 2 D . In order to solve the pole equation ( |1.2| ) to this 
accuracy, we actually need the electrostatic polarization operator n o(0,p) for p 2 = —m 2 
up to order e 4 . That is, we still have to evaluate the momentum dependence of n o(0,p) 
with the prescribed accuracy. Because of the restriction to soft momenta (p = im ~ eT), 
this can be easily done: At one- loop order, we have 11^(0, p) = — m 2 — a e 2 p 2 + 0(e 2 p 4: /T 2 ), 
with the coefficient a computed in App. A (see eqs. (p.7| )- (|B.8|) ). At two-loop order, only 
the momentum dependence of the static diagrams is relevant for Iloo(0,p) to order e 4 (for 
non-static graphs, the leading non-trivial behaviour at small momenta p -C T is ~ e 4 p 2 , 
which is of order e 6 T 2 for p ~ eT). The only static two loop graph which is momentum 
dependent is the diagram l.c, which gives the following contribution to n 00 (0,p): 

C(p) = 4 m< r / ( dk) / (clq) ;, + m , ^ (k + q) l + M2 - («o> 



(For p = 0, this obviously reduces to Cq , eq. ( |4.46| ).) It is easy to see that the momentum 



dependence of Cq (p) does matter to order e 4 if p ~ eT: even if eq. Q4.60| ) is indeed 
analytic with respect to p 2 as p — * 0, the thermal masses which ensure this analitycity 
are precisely of order eT; it follows that, for momenta p ~ eT, all the terms in the small 
momentum expansion are of the same order and the exact expression for C 2 (p) should 
be used. 
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When the previous results are set together, we obtain, for p ~ eT, 

n o(0,p) = -m 2 -ae 2 p 2 - 5Z 3 p 2 + 11$ (0, 0) + 11$ (0, 0) (4.61) 

+ (C 2 Q (p)-C 2 (0))+O(e 5 T 2 ). 

The counterterm 

iz > = Sri (4 ' 62) 

accounts for the photon wave-function renormalization in the minimal substraction scheme 
and cancels the divergent piece of a in the limit D — > 4 (see eq. (p.8| )). This leaves the 
following renormalised value for the coefficient a: 

<ls =3(^( 7 + 2 + " 1 4^) (4 ' 63 > 
To order e 4 , the solution to eq. (|1.2|) is 

m% = m 2 (l - a R e 2 ) - 11$ (0,0) - 11$ (0,0) - (C 2 (im) - C 2 (0)) + 0(e 5 T 2 ). (4.64) 
The difference Cq (p = im) — Cg (0) is readily computed from eqs.( |4.6(J( ) and (|4.46|) as 

VT\ 2 r°°dx _ (2M+m)x ( smhmx 



CL(V - *») - C(0) ={^) J„f ^— -1 - 1 J (4.65) 

(e 2 T\ 2 ( M 2M + m M + m, 2(M + m)\ 
In — — — In ■ 



\ 2tt / y m 2M m 2M + m J 

By also using eqs. (|4.9| ) and (|4.59|) for 11$ (0,0) and 11$ (0,0), we finally obtain 

3 2 a 3 a A a 4 



m 2 =T 2 [ e --^-+b- 1 



D \3 4vr (2tt) 2 (12tt) 2 

where b is the positive number 



+ 0(e 5 T 2 ), (4.66) 



b = - 1 + 57S + { l + t) ln t 1 + ^) = 07211328 " • (467) 

and e is the gauge coupling at the scale p (e = e(/i)), as defined by the minimal substrac- 
tion scheme. The expression (|4.66|) is independent of the choice of p if de/dln/i = e 3 /A8ir 2 , 
which is precisely the one-loop (3— function of scalar QED. To verify that the screening 
mass (|4.66|) is actually gauge independent, we rely on the gauge invariance of n o(0,0) 
and on the fact that the momentum dependent terms which enter Hoo(0,p) to order e 4 
(see eq. ( |4.61| )) arise from gauge- independent diagrams. (The latter assertion is obvious 
for the one-loop contribution, I1qq(0,p); as for the two-loop contribution of eq. (|4.65| ). 
this only involves the time-time piece of the static photon propagator, which is the same 
in all covariant gauges.) 



26 



4.6 Summary: effective theory for static modes 



In the previous sections, we have computed the zero-momentum limit (u = 0, p —>■ 0) 
of the photon polarization tensor to order e 4 . There are two type of terms occuring in 
the perturbative expansion, and we want to discuss this in more detail now. Consider 
the electric piece, n 00 (0,0), as given to order e 4 by eqs. (|2.5|), ( |4.9|) and ( |4.59| ). Some 



of the contributions in these equations arise from purely non-static one- and two-loop 
graphs, which are computed with the standard Feynman rules (i.e. without resummation), 
except for the fact that the terms with u n = are excluded from the sums over Matsubara 
frequencies. This is the case for the leading order electric mass, ( |2.5|) , and also for the two- 
loop contribution of eq. ( |4.57| ) . Since no IR problem arises in such diagrams, naive power 
counting applies. Quite generally then, n-loop diagrams with no static internal line will 
contribute to n o(0, 0) to order e 2n . The second type of terms involve one- and two-loop 
diagrams in which all the internal lines are static and include the corresponding thermal 
masses. The non-analytic contribution of order e 3 , eq. ( |4.9|) , belongs to this category, but 
this is also the case for the order-e 4 contribution of eq. ([4.58|) . As already mentioned, 
the mixed two-loop graphs (where one of the internal frequencies is vanishing, and the 
other is not) contribute to order e 5 . In general, it can be verified by power counting that a 
n-loop diagram (n > 1) with only static (dressed) internal lines contributes to n o(0, 0) to 
the order e n+2 , and eventually dominates over the corresponding non-static contribution 
(of order e 2n ) as soon as n > 3. As for the mixed n-loop graphs (with n > 2), their 
leading contribution to n o(0, 0) is also of the order e n+2 for n > 3 (it is only of the order 
e 5 = e n+3 for n = 2). 

A different, probably more transparent, way to look at the perturbative expansion 
after the resummation of the thermal masses is to consider, as an intermediate step, the 
effective three-dimensional theory for static and long wavelength (p ^ eT) fields. This is 
obtained after integrating non-static loops with static external linesp0[- The lagrangian 
of the effective theory reads (_Dj = Qj — ie 3 Ai) 

£eff = ifg + ^(<^) 2 + ~ + ~ ml Al + (4-68) 

+(A0) t (A0) + M o 2 0t0 + SC. 

The magnetostatic gauge fields A(x), i — 1, 2, 3, the electrostatic (gauge-invariant) field 
A (x), and the scalar field 0(x), may be identified, up to normalizations, with the zero- 
frequency components of the original fields. The term 5C contains higher order vertices, 
but also derivative corrections to the n-point vertices shown explicitly; in particular, it 
contains the counterterms necessary for UV renormalization. All these operators, as well 
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as the parameters and Mq, are obtained as power series in e 2 by evaluating non- 

static diagrams with static external lines in the original theory, and by expanding with 
respect to the external momenta. To leading order, e| ~ e 2 T, m 2 , ~ m 2 and Mq « M 2 . 
Generally, a vertex with n static external lines which is absent at the tree-level is first 
induced at the order e n , via non-static one- loop graphs. Actually, some particular vertices 
may be induced at a level higher than e n , or may even vanish, because of some symmetry. 
The lagrangian (|4.68|) is invariant under static gauge transformations. 

The loop corrections in the effective theory generate all the static diagrams of the 
original theory. They also include the mixed graphs: indeed, a non-static subgraph of 
an original mixed diagram appears as a bare vertex in the effective theory, while a static 
subgraph appears as a loop correction. The calculation that we have done before can be 
understood simply in terms of this reorganized perturbation theory. Consider for exemple 
the calculation of n o(0, 0). In the effective theory, we have n o(0, 0) = — m 2 , — XU (p = 0), 
v where E^ (p) is the self-energy of the field A in the effective theory. Thus, to get 
noo(0,0) to order e 4 , we need both m§ and £U (0) to order e 4 . We already know the 
result for m 2 ,: by evaluating the corresponding non-static one- and two-loop graphs, we 
have obtained (see eqs. Q2.5Q and (|4.57|) ) 

-0 = ^{l + ^e 2 +O(e 4 )}. (4.69) 

To order e 4 , XU o (0) is given by the one- and two-loop diagrams in Fig. 8, where the 
leading order vertices (e| = e 2 T) and massive propagators should be used (Mq = M 2 and 
rriQ = m 2 ). Their evaluation is straightforward, and leads to 





S A (0) = -—-+ _ l + 41n 1 + -= , (4.70) 



which reproduces the contribution of the relevant static graphs in the original theory. By 
adding the contributions ( [4.69j ) and ( |4.70| ), we recover our previous result for n o(0,0). 

It is easy to see the systematics of the higher order corrections. The perturbative 
expansion of il 00 (0,0) decomposes into the sum of two terms: the first one, m 2 ,, is a 
power series in e 2 and arises from the non-static diagrams of the original theory; the 
second term, ^a {p = 0), is obtained as a power series in e in the effective theory. The 
fact that it is e rather than e 2 which governs the perturbative expansion in the effective 
theory is due to the fact that the loop integrals in the effective theory would be IR 
divergent in the absence of the thermal masses. Thus, a n-loop diagram contributing to 
Y,a has an explicit factor of e| n ~ e 2n T n ; since it has dimension two, it is proportional 
to e 2n T 2 (T/M) n - 2 ~ (e 4 T 2 )e n_2 , where M ~ eT is any of the bare masses m or M , 
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or the external momentum, playing the role of an infrared cut-off. Since n > 1, we see 
emerging a power expansion in e, starting at the order e 3 . We have assumed here that all 
the vertices entering the n-loop diagram are of the type explicitly depicted in eq. ( |4.68| ), 
i.e. vertices which are already present in the original theory at the tree-level. Actually, 
the new vertices contained in 5C will also contribute, but only starting at the order e 5 (as 
may be verified by power counting). In particular, to obtain the electric screening mass 
to order e 5 we need the bare parameters Wq and Mq to two-loop order in the original 
theory (e.g. eq. (|4.69| ) for wig) and the corrections to S J 4 (p = 0) up to three-loop order 
in the effective theory. The latter will involve one-loop diagrams with some of the 4-point 
effective vertices from 5C. 



5 Conclusion 

In this paper we computed directly the Debye screening masses for QED and SQED two 
orders above their lowest nontrivial values and verified explicitly that the corresponding 
magnetic masses vanished in the same approximation. It was also shown that the results 
were both gauge and renormalization group invariant. For QED, we further argued that 
the vanishing of the magnetic mass holds to all orders in perturbation theory: Since the 
(odd) Matsubara frequencies of the fermionic propagators play the role of a mass, the 
static photon lines are always coupled to massive fields; there is no closed loop of massless 
fields. This is enough to ensure the analyticity with respect to p 2 (p being the external 
three momentum) of the Feynman diagrams contributing to the static polarisation ten- 
sor. These arguments can be extended to SQED as well, in spite of the fact that the 
propagator of the scalar particle may be at zero Matsubara frequency: Once thermal 
masses of order gT are given to the charged particles and to the longitudinal photons by 
the appropriate resummations, the situation in SQED becomes similar to that of QED 
— - the magnetostatic fields, which remain massless, couple only to massive particles. By 
contrast, the properties of the polarisation tensor are different in non-abelian theories. 
For example, in QCD there is no identity like fll.4j) and, in fact, the self-interactions of 
the magnetostatic gluons are believed to generate a magnetic screening mass ~ g 2 T |24[ . 

In QED, the above arguments on the analyticity, with respect to p 2 , of the static 
polarisation tensor are valid for momenta p smaller than T because the odd fermionic 
frequencies ~ T provide the relevant infrared cut-off. This is why, for example, in eq. 
(3.21) we could use the small momentum expansion of IToo(0,p) up to p ~ im. On the 
other hand, in SQED, the infrared cut-off is provided by the screening masses ~ gT, and 
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the small p 2 expansion holds only for p < gT. Thus for example, in the calculation of the 
electric mass at order e 4 , we could not expand II 00 (0,p) for p ~ im (see eq.(4.65)). In 
fact singularities do occur for imaginary values of p of order gT. For exemple, *IL.j(0,p) 
in eq. ( f4.12|) has branch point singularities at p = ±2iM, and *£(0,p), eq. ( }4.14|) , has 
similar singularities for p = ±iM. It is interesting to note, however, that while *Uu(0,p = 
±2iM) = 0, *S(0,p = ±iM) diverges. This latter behaviour prevents the calculation of 
the correction of order e 3 T 2 to the screening mass M 2 of the scalar field along the lines 
followed in this paper for m 2 D . 

An analytical structure similar to that of the scalar self-energy *S(0,p) has been 
observed in the calculation of the non-abelian Debye mass at the next to leading order, i.e. 
2 D up to order g 3 T 2 [[3], j23fl . This analogy between the scalar self-energy in SQED and 



the gluon self-energy in QCD can be expected from the similarity between the charged 
scalar sector of £ e //, eq. ( |4.68| ), and the electrostatic sector of the corresponding three- 
dimensional effective theory for high temperature QCDpO], |21fl . In Refs. 0, |23]| , the 
logarithmic divergence in the electrostatic gluon polarisation tensor at p = im has been 
cut-off by giving the static gluon a magnetic mass of order g 2 T. Because of the analogous 
difficulty in SQED where no magnetic mass can be generated, we feel that it is worthwhile 



to look for alternative treatement of this problem^ 



We had stated that the derivation of ( |1.3| ) was formal but checked its correctness 
explicitly to fifth order in QED. On the other hand if one has faith in that identity, then 
our verification may be reinterpreted as actually having provided a nontrivial check on 
the correctness of our resummed perturbation expansion. For SQED we did not verify 
the identity to the fourth order but in this case we repeated our computations using two 
different resummation schemes, thus providing again some useful cross checks. 

Of course, static correlators are not the only relevant probes of a plasma. The 
plasma frequency and damping rate are two important quantities which can be deduced 
only from dynamic correlation functions. Compared to static quantities, the computation 
of dynamic quantities is more involved — not only must one work in real-time (that is an 
analytical continuation of imaginary time, or the real-time formalism itself), but one also 
has to use the full machinery developed in 0] in the resummation necessary for higher 
order calculations. As far as we know, of the four possible explicit resummation methods 
mentioned in the Introduction, only method (c) (use of rearranged lagrangians incoporat- 
ing the hard thermal loop effects) has so far been used for dynamical calculations beyond 
leading order 0, [T7|, E3, [III. 
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Appendix A 

In this Appendix, we review a derivation of the hard thermal loop for the photon 
polarization tensor which emphasizes its classical character and the ensuing universality 
of its structure for any gauge theory. To start with, we recall that IT^ = 5j l ™ d /5A U ', where 
jind jg curren t induced by the electromagnetic field. To leading order in e 

r c] 4 K r d 3 k 

jf d {x) = e J —2K>W{K,x) = e J - [6N + {k,x) - 5N_(k,x)} , (A.l) 

where W(K,x) = 2ir 5(k^ — el) [8(k )5N + (k, x) + 9(— k )5N_(— k, x)} is a gauge-invariant 
Wigner function for the charged particles, ultimately related to the scalar propagator in 
the presence of the gauge fields fllj. The 5-function defines the mass-shell for the scalar 
particles; in the present approximation, it coincides with the free mass-shell, = el (with 
eu = k for massless particles). The equation satisfied by the Wigner functions 5N±(k,x) 
is 
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(v ■ 8 X ) 5N ± (k, x) = T e v • E(x) (A.2) 

de k 

where = (1, v), v = k/e^, v ■ d x = d t + v • V, E is the mean electric field, and 
No(ek) = l/(exp(/3efc) — 1) is the equilibrium occupation factor for bosons. Note that for 
an abelian plasma this is merely the (linearized) Vlasov equation. By combining eqs. (A.l) 
and (TATg ), one finds jf d (uj, p) = U^u, p) A"(lu, p), with 

n, v (u, p) = 2e J — — [6» 5„ - v,v v —— j . (A.3) 

For e,t = k, the radial integral over k can be easily performed, and one recovers the 
expression (|2.4|) . 

In abelian plasmas, j^ nd is gauge-invariant and linear in A^, as shown by eqs. ( A.l|) 
and ( |A.2| ), so that there is no hard thermal loop for the ra-photon vertices with n > 2. 
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In contrast, the QCD induced current is a color vector, i.e. transforms in the adjoint 
representation of SU(N). Then, as required by gauge symmetry, the kinetic equation 
(|A.2|) involves a covariant line- derivative (i.e. v ■ d x is replaced by v ■ D x in its l.h.s.), so 



that the color current is a non-linear functional of the gauge potentials. It follows then 
that the higher functional derivatives of j" with respect to A a ^ are also non- vanishing, and 
define hard thermal loops for vertex functions with an arbitrary number of soft external 
gluon lines. However, for weak fields, \A P \ — >• 0, the kinetic equation ( |A.2|) is formally the 



same for abelian and non-abelian plasmas. This explains why the expressions obtained 
for the polarization tensor in this approximation are, up to a trivial factor which counts 
the relevant degrees of freedom, identical for QED, QCD or scalar QED. 

As for the hard thermal loop for the scalar self-energy, this describes the response 
of the plasma to long-wavelength scalar mean fields. Since according to eq. (O), E^ 2 -* is 



a local operator, gauge covariance requires this operator to be independent of the gauge 



mean field |TT|], so that there are no hard thermal loops for vertices involving scalar and 
photon external lines. We conclude that the self-energy corrections in eqs. Q2.4Q and 
Q4.2D are the only hard thermal loops for SQED, in accordance with the power counting 



arguments of Refs. H, 10 



The previous derivation of IL^ can be generalized to take into account the thermal 
mass M aquired by the charged particles. Since M <T, the only effect is the modification 
of the mass shell condition which becomes k$ = e 2 = k 2 + M 2 . After inserting this into 
eq. (|A.3|), one obtains 



IV(cu = 0, p) = -m 2 (T, M) <V)<W (A.4) 



with 



2,™ 2 r d3k dN o e 2 T 2 e 2 MT ~/ 2%r2 \ ^ r^ 

M > = - 2e 7 w? -d = ~3 sr + ° ( ) • < A ' 5 ) 

As advertised in section 3.2, this is the correct electric mass up to the order e 3 . 



Appendix B 

In this appendix, we derive some formulae which are referred to in Sect. 3. We 
consider first the small momentum behaviour of the one-loop static polarization tensor, 
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n»^(0,p). For the magnetic sector, the diagrams in Fig. 3 imply 



n^(0,p) = (e^) 2 /[dK] 



(2k + p) 



2p-i) 

u; 2 + k 2 ( W a + k 2 ) (u, 2 + (k + p)2) 



(b.i; 



The dominant contribution as p — > is given by the u n = term in the Matsubara sum 
and it is linear in p: 



n^(o, P ) 



(e//) 2 P 2 T/(dk)- 



n=0 



k 2 (k + p) 2 



(B.2) 



e 2 pT r°° dx 



An 2 Jo x 
For the time-time component we obtain 

2 



In 



x + 1 



x — 1 



e 2 pT. 



I4 L (0, P ) = -{e^ f J[dK] 



Au 2 



c 2 + k 2 (a ; 2 + k 2 )(u; 2 + (k + p) 2 ) 



(B-3) 



-m 2 + {e^f J [dK\ 



Awl 



1 



;^ 2 + k 2 ) v^ 2 + (k + P ) 2 ^ + k 2 , 



where 



m 



-njf (0, 0) = 2(D - 2) (e//) 2 | [dQ] Sb(Q) 



e 2 T 2 



(B.4) 



Note that the static mode uj n = does not contribute to the last sum in eq. ( |B.3j ); accord- 
ingly, one can expand the denominator for small p without generating IR singularities: 



1 



1 L P 2 + 2k-p (p 2 + 2k-p) 



+ k 2 



(^ 2 + k 2 ) 



u; 2 + (k + p) 2 oo 2 n + k 2 
By keeping only the terms quadratic in p, we obtain 

II#(0,p) = -m 2 -ae 2 p 2 + 0(eY/T 2 ), 



(B.5) 



where 



A^ 2e J [dK] 
5-D 



to. 



[uj 2 + k 



2\3 



D-lu 2 + k 2 



vB.6) 



(B.7) 



Ai* [dK]'S 2 (K). 



and some integrations by parts have been necessary to get the last line. By also using 
eq. (|3.7| ) with n = 2 and by expanding in e, we obtain 

,2 \ 



1 



3(4tt) 2 



;+7 + 2 + ln^ + 0( e ) 



(B.8) 
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Finally, we turn to the evaluation of the integrals ( |4.32| ) and (|4.46 ) which contribute 



to n 00 (0) to order e 4 . We organize this computation as follows: 

A 3 + C oo = a 1 + a 2 + a 2 , (B.9) 

where 

a, S 2 W T 2 /(dk) /(dq) j^l^ + 7F ± Jp - , (B.10) 

a. 2 , 8 M ' M 2 T 2 /(dk, /(dc) i (j^p^ (q + k) l + M2 . (B.11) 



and 

a 3 = -4 (e//) 4 T 2 / (dk) / (dq) - 1 - ro tttt f 4 ~ o 1 J • (B.12) 

d VP; J K 'J y H; k 2 + M 2 (q + k) 2 + M 2 \q 2 q 2 + m 2 / v ; 

All the integrals above are UV and IR finite in D = 4, and thus we can set e = 0. The 
evaluation of ai is straightforward, with the result 

fe 2 T\ 2 m + M , 

In order to compute a 2 and a 3 , it is convenient to use the coordinate-space representation 
of the static propagators, by writing 

1 r - P ~ Mx 1 1 /• 

- = / d 3 xe ik S - - = — — / d 3 xe ik S e~ Mx (B U) 

fc2 + M 2 J axe 47rx ' (P + M 2 ) 2 MM J ' 1 j 

with x = Then one obtains 

(e 2 T\ 2 

a 2 = 2^— j , (B.15) 

and 

- ^ " 4 (£ ) 2 f f (<~ " ^ - (£ )" - ^ < B -> 

By adding together eqs. ( p.!3|) , ( |B.15|) and (|B.16|) , we obtain the expression (|4.58|) . 
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Figure captions 



Figure 1. The skeleton diagrams for the self-energy of the photon in scalar QED. 

Figure 2. The two loop contributions to the photon self energy in QED. 

Figure 3. Two loop diagrams in QED in which the photon internal line is dressed with 
the thermal mass. 

Figure 4. One loop contributions to the polarisation tensor in SQED. 

Figure 5. One loop contributions to the scalar self-energy in SQED. 

Figure 6. One loop contributions to the polarisation tensor in SQED, with the static 
internal lines dressed by thermal masses. 

Figure 7. One loop contributions to the scalar self-energy in SQED, with the static 
internal lines dressed by thermal masses. 

Figure 8. The self-energy of the static photon in the three-dimensional efffective theory 
(see eq.(4.70)). Diagram (c) gives no contribution in dimensional regularisation. Dashed 
line: longitudinal photon. Wavy line: transverse photon. Full line: scalar field. 
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